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Abstract 



■ A general method for reducing tensor form factors, that appear 

in one-loop calculations in dimensional regularization, to scalar in- 
tegrals is presented. The method is an extension of the reduction 
scheme introduced by Passarino and Veltman || and is applicable in 
all regions of parameter space including those where kinematic Gram 
determinant vanishes. New relations between the the form factors that 
valid for vanishing Gram determinant play a key role in the extended 
scheme. 



1 Introduction 

The evaluation of one-loop tensor integrals, that was pioneered by Brown 
and Feynman || and systematized in the work of Passarino and Veltman |IJ , 
is by now a mature science. The latter authors defined a set of generic tensor 
form factors in terms of which any one-loop tensor integral could be written. 
Thus, for example, the 3-point tensor integral with two Lorentz indices, 



is written in terms of generic tensor form factors Cij{p\ 1 p\ 1 p\\m\ 1 m 2 , 1 m\) 
defined by 



where p$ = Pi+p 2 and we have omitted the arguments of the Cy's in eq.(0). 
This choice for defining the Cy's has the obvious advantage of its great 
simplicity. Passarino and Veltman showed how the C^-'s and their analogous 
4-point form factors, Dij could be reduced to scalar integrals and 2-point 
form factors. 

In ref. it was shown how to reduce the Cy, to unique form consist- 
ing entirely of scalar integrals. The advantages of doing this are manifold. 
The resulting expressions are often simpler in form, containing just a few 
scalar integrals with coefficients that are rational functions of the masses 
and momenta. All non-trivial analytic structure, such as cuts in the complex 
plane, are isolated in the scalar integrals which can be evaluated separately 
taking care of these technicalities 0, [|. Having a general and well-defined 
procedure for reducing tensor form factors to scalar integrals lends itself to 
automation and means that quite lengthy calculations can be undertaken. In 
such calculations it is essential that stringent checks be performed as to the 
correctness of the final results. These require the comparison of two algebraic 
expressions whose equality can only be determined if both have been fully 
reduced to a unique form containing only scalar integrals with no remaining 
tensor form factors || |7j . 




in 2 [q 2 + ml] [(q + pi) 2 + m§] [(q + pi + p 2 ) 2 + m|] ' 



(1) 



C flu (p 1 ,p 2 ;m 2 1 ,ml,ml) 



Pl M Pl M C 2 i + p 2 ^p 2u C 22 + {VXixViv + PluP2^}C 23 + S^C 2 4 (2) 
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Passarino and Veltman applied their reduction scheme numerically to the 
to the process e + e" — > fi + yT and in doing so highlighted a significant dis- 
advantage of their approach. It is sometimes the case that strong numerical 
cancellations occur during reduction. These authors, for example, found it 
necessary to calculate to 40 significant figures in order to obtain sufficient 
accuracy. Performing the reduction algebraically as in ref . || has the advan- 
tage of eliminating many of the strong numerical cancellations and hence 
leads to improved numerical stability. It can also lead to rather compact and 
convenient expressions for physical observables. 

Still the reduction scheme of Passarino and Veltman has other vexing 
problems. The reduction formulas for 3- and 4-point functions involve the 
inversion of kinematic matrices and therefore it breaks down when the deter- 
minants of these matrices, Gram determinants, vanish. In ref. @ it was shown 
how to extend the Passarino- Veltman reduction scheme to treat situations 
of vanishing Gram determinants. This was based on the observation || that 
in such cases, the n-point form factors reduce to (n — l)-point form factors. 
This extended reduction scheme was implemented in the REDUCE program 
LERG-I ||. Even here, the procedure broke down for certain special cases 
of external momenta and internal masses. In ref.|| extensions were made 
to enlarge the region of parameter space for which the reduction could be 
performed but some problematic holes nevertheless remain. The extended 
scheme was implemented in both REDUCE || and Mathematica [ 1Q"| . In 



this paper we present an alternative reduction scheme for the case of vanish- 
ing Gram determinant derived directly from the original Passarino- Veltman 
approach. New techniques involving differentiation with respect to mass ar- 
guments are introduced that make this and earlier schemes viable everywhere 
in parameter space, except perhaps where true infrared divergences or mass 
singularities appear. 

Although expressions that have been algebraically reduced to scalar in- 
tegrals exhibit improved numerical stability over those in which the tensor 
form factors are calculated numerically in a hierarchical fashion, they can still 
become problematic as kinematic determinants become small. This problem 
was addressed by van Oldenborgh and Vermaseren [I t] who employed a dif- 
ferent tensor basis, constructed from so-called inverse vectors, to define the 
form factors. This approach has the effect of simplifying the contraction 
of Lorentz indices and yields expressions that can exhibit improved numeri- 
cal stability for small Gram determinants. Since the form factors that they 
obtain are linear combinations of those defined by Passarino and Veltman, 
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their method can then be viewed as a way of arranging expressions to reduce 
numerical instabilities but it does not show, as they state, that the approach 
of ref . P| "is most suited for numerical application" . 



Ezawa et al. [12] have performed the reduction of the tensor form factors 



using an orthonormal tensor basis. 

Pittau |H| has proposed a reduction method, based on 7-algebra, that 
avoids the appearance of Gram determinants and is applicable for loop inte- 
grals with at least two massless external legs. 

More recently Campbell et al. [|l4j have pointed out that the Gram de- 
terminants are artifacts of the procedure used to reduce tensor form factors 
to scalar integrals. They show how to write Feynman parameter represen- 
tations of tensor integral form factors in terms of scalar integrals in higher 
dimensions and derivatives of ordinary scalar integrals. The various Feynman 
parameter integrals are separately finite and highly stable as the condition 
for vanishing Gram determinant is approached. In practice, however, evalu- 
ating the Feynman parameter integrals yields terms between which there is 
eventually a strong additive cancellation and numerical instabilities do set in 
although their appearance is much delayed. 

Other tensor integral reduction schemes have been put forward. In quan- 
tum gravity, the computation of tensor integrals by reduction to scalar in- 
tegrals was pioneered by Capper and Liebbrandt |]T3j| both for the graviton 
self-energy and the fictitious particle loop. Davydychev has ]16| obtained 
reduction formulas from relations between integrals in different numbers of 
dimensions. Recursion relations obtained recently by Tarasov [|17) have com- 
pleted the picture to make this a practical approach but here again problems 
arise when the Gram determinant vanishes. Bern et al. [|11J derive tensor 



integrals by differentiation of Feynman parameter representations of scalar 
integrals and introduced many of the basic techniques employed in ref. [|T1|] . 

All reduction schemes have certain strengths and weaknesses and may 
be more conveniently applied in certain circumstances than others. In the 
scheme we present here we take the point of view that one important advan- 
tage of Passarino-Veltman reduction is the simplicity in the way the form- 
factors are defined (c.f. eq.(pj)). This simplicity may come at the price of 
numerical instability near phase space boundaries but even there appropri- 
ate Taylor series expansion of the scalar form factors may eliminate it. In 
such cases other methods which employ more numerically stable form factors 
may be more easily applied. The methods given here can also lead to lengthy 
algebraic expressions when there are many different masses and momenta in 
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the problem. These expressions can be costly in terms of computer storage 
and processing time although such issues are of ever-declining importance. 

In this work we consider the form factors, CV,- and Dij, as well-defined an- 
alytic functions of their arguments everywhere in parameter space. Feynman 
parameter representations exist for these functions which, can be expected to 
be analytic, at least in the physical region, except for isolated singularities. 
The reduction scheme is one way of obtaining an expression for, or value of, 
these functions in a given region of parameter space. Note that this means, 
for example, that C21 and C22 with arguments C^p 2 ,^, 0; m^, m^, m^) re- 
main well-defined entities even though, with p± — p 2 , eq.(fj) does not uniquely 
specify them. It also means that relations such as (|2]) are valid everywhere 
in the parameter space. This approach lends itself well to implementations 
with computer algebra. 

Our aim here is to obtain valid analytic expressions for the tensor form 
factors in terms of scalar integrals for all regions of the parameter space 
where those form factors do not experience genuine singularities. We do 
not address the problem of numerical stability, since it has been intensely 
studied elsewhere, but concentrate rather on on the region where the Gram 
determinant exactly vanishes. A new ingredient is the use of additional 
relations between tensor integrals and scalar integrals that are valid only for 
vanishing Gram determinant. These relations actually follow directly and 
straightforwardly from the original reduction formulas given by Passarino 
and Veltman. 

At first sight it might seem somewhat strange that there is such difficulty 
in evaluating the tensor form factors for vanishing Gram determinant. After 
all one could, in principle, vary one of the arguments away from region of 
vanishing Gram determinant and then use l'Hopital's rule to take the limit. 
One finds, however, that this procedure regenerates a vanishing Gram deter- 
minant in the denominator of the resulting expressions. The new relations 
between the form factors allow l'Hopital's rule to be applied should it be 
necessary in extremely singular situations. 

2 Notation 

In dimensional regularization the dimensionality of space-time, n, is contin- 
ued away from n = 4. Quartic, quadratic and logarithmic divergences corre- 
spond to poles at even integer n > 0, 2 and 4 respectively. One introduces 
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the logarithmically divergent constant, 



n _ / n 



A = ^Y(2--) (3) 

-2 . , 

— 7 — In 7r (4) 



n- 4 

as n -> 4 where 7 is Euler's constant. 

Note that it is essential to keep careful track of background terms that 
may be generated when A is multiplied by polynomials in n. For example 

(n - a) A -»• -2 + (4 - a) A 

as n — > 4. 

A detailed exposition of dimensional regularization is to be found in 
refs pi [TP| along with results for certain standard integrals. For our pur- 
poses we need in addition, 

n . 

n PkPxPhPv »tt? 1 ) v ( _ 1 }\i c i c i c i c 

P (p* + 2k-p + m 2 )<* " (m 2 -P) Q -ir(a)\ v" 2J K A » v 

+ Ir(«-l-|)(m 2 -A; 2 ){^} K A^ (5) 
+ ir(a-2-|)(m 2 -P) 2 W K ^} 

In defining tensor integrals we shall adopt the notation of ref . |T| in which 
braces denote symmetrization with respect to Lorentz indices. Thus for ex- 
ample, 



{PS} fj.ua = Pn$iya+PJna+Pa$nv (6) 
{b$}fA,uaP = 5[i.u5af3 + S^ a S u /3 + 5^p5 ua (7) 

We will use the notation B (i,j) and Co(i,j, k) in context to denote the 
two- and three-point scalar integrals constructed from the z-th, j-th and fc-th 
denominators of some tensor integral. 
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3 One- and two-point form factors 

Following Passarino and Veltman [|IJ we use a metric in which the squares of 
time-like momenta are negative. 

The one-point scalar integral, A(m 2 ), is quadratically divergent and de- 
fined as, 



A (m 2 ) = /?f— (8) 



71 _o 

7r^ _ / n 



V 2 



«r i-77 ( 9 ) 



The two-point tensor form factors, Bij(p 2 ; m\ , m 2 ,), are defined in Ap- 
pendix A. These form factors are not independent of one another and may 
all be reduced to the form, 

+/?i£o(0; m 2 , ml) + (3 2 B (0; m 2 , m 2 ) + (3 (10) 

in which a.\ and the /3's are rational functions of the momenta and masses, 
p 2 , m 2 , ml, 

B (p 2 ;ml,ml) = [ , f 2TT7^ — 7^ 27 ( n ) 

and 

i?o(0; m 2 , m 2 ) = A - In m 2 (12) 

The coefficient ot\ is non-zero only for quadratically divergent Bij. The 
form factor B43 is quartically divergent and requires the introduction of an 
additional term but we do not make use of it explicitly in the the present 
work. In addition, the partial derivatives of any form factor B^ with respect 
to any of its three arguments can be written in the form (|TT)|) . 

By performing the usual Feynman parameterization of the integral (|TT ) 
one obtains, 

B a (p 2 ;m 2 ,ml) = (-)V?- 2 r(2 - -) / dz z l \az 2 + bz + c]?~ 2 (13) 

^ Jo 

= (-)* f dz /[A-ln^az 2 + bz + c)] (14) 
Jo 
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where B m = B , 5 n = B\ and a = —p 2 , b = p 2 — ml + mf, c = — ze. 
Eq.([14]) follows from eg. (p~3|) by dropping terms of 0(n — 4) and higher. Such 
terms will be discarded without further comment in what follows. In the 
region of n = 4 these integrals become 

Wi ro? , mi ) = 1z)1{a - Ko + 6 + c) + [ } (15) 

The integrand can be decomposed as, 

z i+1 (2az + b) A (ft* + A,)(2as + &) 

: = > OtjZ J H — (16) 

az 2 + bz + c ^— ' az 2 + te + c 

Carrying out the integration for all but the term proportional to (3% and 
comparing with (IT) and (|T5D for i = yields, 



Bn{p 2 ;ml,ml) = 7-7^ Pi B o(p 2 ] rn%, ml) + f3 B (0; m\, m\) 



i 1 1 

+(1 - A) - /3i)-B Q (0, "is. "^) + ~Ti aj I (17) 



3=0 

Full expressions for the form factors B^ can be found in ref. || and will not 
be repeated here. 

The partial derivative of the Bn(p 2 ; ml, m|) with respect to p 2 form fac- 
tors arise in practice in wavefunction renormalization factors and those with 
respect to the mass arguments are needed in intermediate steps in the reduc- 
tion to scalar integrals at kinematic boundaries. Because of the reduction 
(|10D , it follows that one only needs expressions for the partial derivatives of 
Bq and it turns out that the partial derivatives themselves take the form 
(pTDj) . From eq.flTo]) one has for the region of n = 4, 

dBn (p 2 ;ml,m 2 2 ) = - f dz f~ Z ] (18) 



dp 2 Jq az 2 + bz + c 

9B / 2 2 2\ f 1 j (1 _ z ) /m\ 

-{p 2 ;ml,m 2 2 ) = - / dz " ■ - ' (19) 



dml ' ' J az 2 + bz + c 

The partial derivative with respect to m\ easily follows since B is symmetric 
in its mass arguments. 
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Writing 



z(z-l) z{2az + b) (2az + b) 

-« +Pi — 5— t — ; hp — 5 ■ - — ; — (20) 



az 2 + bz + c az 2 + oz + c az 2 + bz + c 



111 



18) for which, 



b 2a + b lb 2 -2ac + ab c 2a + 6 

a o = — T9 ^ — > Pi = i > Po 



ab 2 — 4ac' a 6 2 — 4ac a 6 2 — 4ac 

and 

1-z z(2az + o) (2az + 6) 

a o + Pi — 5—1 — ; ^Po — 5—1 — ; — (21) 



az 2 + bz + c az 2 + bz + c az 2 + bz + c 



in fll9|) for which, 



2(2a + 6) a b + 2c _ 2a + 6 

a o = "T5 ; — , Pi = -75 : — > Po ~ 



6 2 — 4ac 6 2 — 4ac 6 2 — 4ac 

the integrals expressions for the partial derivatives are easily reduced to the 
form (|T0[). Full expressions for a 0) A, Po an d — (p\ + po) in terms of p 2 , m 2 , 
m\ are to be found in Appendix A. 



4 Three-point form factors 

The 3-point scalar form factor Co(pf, p\, p\\ m?,m|, m 2 ) is a function of three 
external momenta squared and three internal masses. Following ref.0 we 
define, 



^ / 2 2 2 2 2 2\ 

d n q 



in 2 [q 2 + mf] [(q + p{) 2 + m 2 ] [(5 + p 1 + P2) 2 + w|] 

with pn = pi + P2- 

This may be written as the Feynman parameter integral, 

n I 2 2 2 2 2 2\ 

2 dzi dz 2 j dz s 5(l - Zi - z 2 — z 3 ) 
Jo Jo Jo 

f d n q 1 

J in 2 [q 2 + pjz 1 z 2 + p\z\Zi + p 2 ZiZ 3 + m\z x + m\z 2 + mlz 3 ] 



(22) 



(23) 
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that clearly displays the symmetry properties of the arguments and from 
which it may be shown that, 



n I 2 2 2 2 2 



, m 

/ * / 2 2 2 2 2 

C {PnP 5 iP 2 ;m 2 ,m 1 ,m 3j . 

„2 „2 „2. _2 _2 



1), 
2\ 



r~i / 2 2 2 2 2 
^Q{P^P^P^ m ^ m ^ m 



2 ) 
2!- 



ri I 2 2 2 2 2 2\ 

C {p 5 iP 2 ,p 1 ;m 1 ,m 3 ,m 2 ) 

ri I 2 2 2 2 2 2> 

Co(pi,pl,pi;m. 



(24) 



are all equal. 

In a similar way as for the two-point form factors, Bij(p 2 , mf, ml 1 ), one can 



define a sequence of three-point tensor form factors, Cy (p 2 ,p 2 ,p 2 



m\,m 2 , 



mi) 



as shown in ref . U . Passarino and Veltman [[y] have shown how to reduce to 
expressions involving Bij and lower Cy form factors. For certain of the 
tensor three-point form factors, this method yields two apparently distinct 
expressions. In ref. [Q] , these expressions were evaluated numerically and their 
agreement was used as a powerful check of the approach. The method that 
they employ can be made use of to algebraically reduce tensor three-point 
form factors to the two- and three-point scalar integrals. 



^-y/222 2 2 2\ 

Ci j {p 1 ,p 2 ,P 5 ;m 1: m 2 ,m 3 ) 



ri l 2 2 2 2 2 2\ 

^a 1 C {p 1 ,p 2 ,p 5 ;m 1 ,m 2 ,m 3 ) 
+PuB (l, 2) + /3 13 B {1, 3) + P 23 B (2, 3) 
+/? 1 B (0;m?X)+/? 2J B (0;m: 



2- 



m Q 



+(3 3 B (0;mlm 2 3 )+(3 



(25) 



where a.\ and the /3's are again rational functions of the masses and momenta. 
The method of ref. [El] involves inverting the matrix, 



X 



pi 

Pl -P2 



Pi ■ P2 

pi 



(26) 



and consequently breaks down with the vanishing of the determinant at 
kinematic boundaries. We will now specify the properties of the external 
momenta, p±, p 2 and p§ for which the determinant vanishes, that is when 

pIpI - (Pi ■ P2? = --M +Pt + Pt~ Ivlvl - 1p\p\ - 2plpl) = (27) 
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If pi = (pijiEi), p 2 = (p 2 ,iE 2 ) and 9 is the angle between their 3-vector 
components then, 



cos9=l w)(w 



± 



1 



El 
\pi\ 



i - 



E 2 
\Pz\ 



(28) 
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From the requirement that cos^ be real and | cos#| < 1 it may be deduced 
that, 

• If one of the momenta p 1 , p 2 , p$ is time-like then the other two are as 
well. Moreover pi = ap 2 for some real constant a. 

• If two of the momenta are light-like then the third is as well. Moreover 
Pi = ap 2 for some real constant a. 

• If one of the momenta is space-like then the other two are either both 
space-like or one is space-like and the other is light-like. In this case 
Pi and p 2 can be strictly linearly independent. 



5 The Derivative of Co 

In the reduction scheme presented in this paper, it is sometimes necessary to 
obtain the expression for the derivative of Cq with respect to one of its mass 
arguments. For completeness the expressions for the derivative of Co with 
respect to any of its six arguments are derived in this section. 

The scalar integral C has a Feynman parameter representation given by 



pi PX 

Co(pi,pl,pi,ml,ml,m 2 3 ) = dx dy 

Jo Jo 



1 



ax 2 + by 2 + cxy + dx + ey + f — it 

(29) 



where, 

a = -j%, b = -pj, c = pl+pl-pl, 

d = m 2 - m 2 +pl e = m 2_ m 2 +p 2_ p 2 ; j = m 2 . 

The derivative of Co with respect to any of it's six arguments is, 

c > = f 1 dx r d ax2x2 + a v 2 y 2 + a *v x y + a * x + a yy + ai ( 30 \ 

J J (ax 2 + by 2 + cxy + dx + ey + f - ie) 2 
where the coefficients, a, are given by 
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Derivative w.r.t. 


a x 2 


a y 2 






(Xy 




pi 





1 


-1 











Pi 


1 





-1 


-1 


1 





Pi 








1 





-1 





2 

m\ 














-1 





m\ 











-1 


1 





m\ 











1 





-1 



It is convenient to define as basis set integrals, Jj, for the decomposition of 
Cq as follows 

f" x(2by + cx + e) , f 1 . f. 2ax + cy + d 
h s ^ ^ rf, _ U s ^ dx j^ dy _ 

, f 1 f 2by + cx + e f 1 f dx + ey + 2/ 
J, = ^ dx j^ iy d2 h = ^ dx^ dy - 2 

where D = ax 2 + by 2 + cxy + dx + ey + / — ie. It is easy to show that the 
Zj's may be written in terms of Co and the derivatives of -Bo's. 



dB (l,3) <9B (1,3) 

dm 2 dm 2 

dB Q {2,3) dB (2,3) 

dm\ dm 2 



The derivative of Bq may be expressed in terms of -Bo's as shown in eq.( |A.7|) . 
Note however that this equation will yield an ill-defined -Bo(0; 0, 0) if m\ is set 
equal to zero. Such terms will generally cancel without any special treatment 



h 
h 
h 
h 
h 
h 



C + 



<9B (1,3) 



950(1,2) 9S (1,3) 



dm 2 

dm\ 
dB (l,2) 

dm 2 
dB (l,3) 

dm\ 

dB (l,2) 
dm 2 



+ 



+ 



dm 2 
dB {2,3) 

dm\ 
dB (l,2) 

dm\ 
<9B?o(l,3) 

dm 2 

dB (l,2) 
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in expressions for physical quantities. It may, in cases where infrared diver- 
gences or mass singularities are present, be convenient to carry a non-zero 
m\ through intermediate steps and set it to zero only in the final results. 
With the above relations in hand the reduction of Cq is straightforward. 

C = fah + f3 2 I 2 + ft/ 3 + Pth + Psh + fol 6 (31) 

The equation for the ft's is obtained by equating the coefficients of x 2 , y 2 , 
etc., in the above. 



/ o 





c 








\ 




( Pi 


\ 


/ 


a x 2 


\ 


2b 


c 


























c 


2a 


2b 













Ps 
















e 


2a 


c 


d 




Pa 










e 


d 





c 


2b 


e 




ft 






a y 




V o 








d 


e 


V J 




V ft 


/ 


V 


Oli 


J 



(32) 



This matrix equation may be solved by first determining (3i, (3 2 , and ft from 
the relation 

c \ / ft \ / a x 2 \ 
2b c ft = (33) 

c 2a 2b J \ (3 3 J \a xy ) 

and then obtaining /3 4 , (3 5 , and ft using, 

2a c d \ ( ft \ ( a x \ ( e \ / ft \ 
c 2b e ft = [ a y - e d ft (34) 

d e 2/ J \ fa ) \ ai J V / V ft / 

It should be noted that the derivative of Cq with respect to any of its 
mass arguments reduces to -Bo's. In this case, a x 2 = a y 2 = oi xy = 0, and 
consequently ft = ft = ft = 0. Since I4, 1$, and Iq can be expressed in 
terms of just -Bo's, it follows that 9C (pf,p|, p\, m 2 , m|, m\)/dm 2 can also be 
expressed in terms of -B 's. 

The expressions for the ft's, for the particular case of the derivative with 
respect to p\, are given below. 

/ ft \ 1 / 2a + c \ 
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M 2bd + cd - 2ae - ce ( f 7 2 °J 



2 



'6 



2(4ab - c 2 )(-bd 2 + cde - ae 2 + Aabf - c 2 /) I cd _ 2 ae 



(36) 



The expressions for the derivatives with respect to the other two momentum 
arguments may be obtained from the symmetry properties of Co (eq.(p4|)). 
The expressions for the for dCo/dm\ are, /3i = fii = (3% = and, 

, / de — 2c/ 

Aaf-d 2 I (37) 




—bd 2 + cde — ae 2 + Aabf — c 2 f \ ^ _ 

Again the derivatives with respect to the other masses may be obtained using 
the symmetry properties of C (eq.(£23|)). 

In case the matrix that needs to be inverted to obtain dC /dm 2 (eq.(|3"4[)), 
has zero determinant, the simple solution is to obtain the expression for the 
derivative with one of the one of the masses shifted away from its value, 
and since the resulting expression will contain only -Bo's, l'Hopital's rule or 
Taylor series expansion may be applied to obtain the limit. Recall that the 
derivative of Co with the most general set of arguments, with respect to any 
of its mass arguments is expressible in terms of -Bo's. The expressions for 
dBo/dm 2 are given in Appendix A. 



6 Reduction of C%j Form Factors 

The Passarino-Veltman scheme for the reduction of all Cij form factors to 
the basic scalar integrals Co and .Bo breaks down when the kinematic or 
Gram determinant, vanishes. The Gram determinant will be denoted, T>, 
and when it vanishes all the CV, form factors, including the scalar integral 
Co, can be reduced to linear combinations of .Bo's as long as they are at most 
logarithmically divergent. Expressions for the quadratically divergent Cy's 
contain the quadratically divergent scalar integral, A. 

The reduction for Co is obtained as follows. The Passarino-Veltman for- 
mula for reducing the scalar integrals Cn and C\i to C 's and -B 's in the 
notation of ref.fll] is, 

(£)-^(£)-K4»Tr)U) ^ 
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where, 



Ri = ^{fiC + B {l,3)-B (2,3)} 
R2 = i{/ 2 C + S (l, 2) -50(1,3)} 

r 2 2 2 

/1 = m l -m 2 - p 1 

e 2 2,22 

f 2 = m 2 -m 3 +p 1 - p 5 

and T> is the determinant of the matrix X, 

V = pIpI-(pi-P2?. (39) 

When V = 0, in order for Cn and C\ 2 to exist, the following equations must 
be satisfied, 

ptRi-pi-p 2 R2 = (40) 
-Pi-PiRi+plR* = (41) 

Substituting for Ri and R 2 one obtains, 

Cq = (pi-p 2 )B (l,2) - (grPg +p|)5q(1,3) +pp?o(2,3) 

r pjB (l,2) - (pl+ Pl -p 2 )B (l,3) + (pi-p 2 )A)(2,3) / 

C ° ~ {PvP2)h-p\f2 (43) 



respectively. Equations (g2j) and (|43|) give explicit expressions for C when 
the kinematic determinant vanishes. Although both of these expressions are 
equivalent, it is useful to keep both expressions at hand because one encoun- 
ters situations where the denominator of one of the expressions vanishes but 
that of the other does not. For example, when p\ = and p\-p 2 = 0, the 
denominator of the expression in fl42|) vanishes but that of (|43"D need not. 

The eq.s ([12]) and ( f!3"D are new formulas valid when D = and represent 
an alternative to the method described in ref.s|3|. [|. 

At this point we introduce some notation to make the equations more 
concise. The cofactor matrix of the matrix X, will be called X. Thus, 

P\ ~Pl -P2 

-Pi -P2 P\ 







X\2 N 


H 




V ^21 


3^22 / 
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We also define, 



Y l = xnh+xvh (44) 
Y 2 = X21/1+X22/2 (45) 



In this notation the two equations for Co, ([Q) and (53), can be written 
concisely as, 

C = ^{-x i2 B (l, 2) + (x l2 - x a )B (l, 3) + x a B (2, 3)} (46) 

where i — 1, 2. These eq.s(S5p run into difficulties when both of the 5^ 
vanish. In what follows it will be shown how to bypass such problems by 
differentiation with respect to mass arguments. It can be demonstrated that 
the vanishing of both Yj's is a necessary condition for the breakdown of the 
method described in ref. || and many of the techniques we introduce here 
can be used to extend that method as well. 

To obtain the reductions for Cn and C\ 2 when V = 0, one could attempt 
to obtain the limit of the rational functions in eq. (|38]) using l'Hopital's rule 
but we find that when the partial derivative of Co with respect to one of its 
momentum arguments, is expressed in terms of Co and the -Bo's, the Gram 
determinant, eq.(^), appears in the denominator and the application of 
l'Hopital's rule in this manner does not yield a limit. Instead, the reductions 
for Cn and C12 when T> = can be obtained by looking at the general 
formulas for the C 2 j form factors, 

° 21 \ = x-i ( Ri \ W lCl1 + 3 ) + Bo & 3 ^ ~ ° 2i 

C 23 J \R 5 J ~ V\l{f2Cu + B 1 (l,2)-B 1 (l,3)} 



R 3 \ 


X 


Rb J 


V 


Ra \ 


X 


Re 


V 



C 2 z \ _ x -x( Ra\ _ Xfl{f 1 C Vi + B 1 {l,3)-B 1 (2,3)} 



C 22 J \R 6 J V\ \{f 2 C l2 - B 1 (l, 3)}-C 



-24 



Hence when T> = 0, we must have, 



xaRa + xuRs = (47) 

xnR4 + Xi 2 R 6 = (48) 

Substituting for R 3 , R4, R 5 , and Rq one obtains, 

YiC n = -x i2 B 1 (l,2) + (x i2 -x il )B 1 (l,3)-x il B (2,3) + 2x il C 24 (49) 

YiC 12 = (x i2 -x il )B 1 (l,3)+x il B 1 (2,3)+2x i2 C 24 (50) 
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Unlike the case for Co where the right hand side did not include any Cy form 
factors, here the logarithmically divergent C 2 4 is encountered on the right 
hand side. The appearance of C24 causes no difficulties because it is given in 
closed form by, 

^24 = \ ~ \{hC n + f 2 C l2 + 2m\C, - B (2, 3)} (51) 

Note that this is not affected by a vanishing denominator when D = 0. 
Substituting this expression for C 24 into (pE9p and (BO) one obtains, 



x i2 h 2Y i + x i2 f 2 \ C 12 J \R 12 



(52) 



where, 



Ru = -2mlx il C -2x i2 B 1 {l,2) + 2{x i2 -x i i)B 1 {l,3) 

-XnB (2,3) +x a 
R 12 = -2m 2 1 x i2 C + 2(x i2 -x il )B 1 (l } 3)+2x il B 1 (2,3) 

+x i2 B (2,3) + x i2 



Once the substitution for Cq is made using eq. ( fig) and B\ is written in terms 
of B [[|, the expressions for Cn and C\ 2 will be given purely in terms of 

Bo's. 

The formulas for the reduction for C 2 j, C$j, and form factors, when 
Z> = 0, are obtained in a similar manner from the general formulas for the 
reduction of the C 3j -, Cy, and C 5 j form factors respectively. The difference 
between the C 2 j and higher form factors, and the ones treated above, C and 
C\j, is that the former have form factors that multiply tensors that have one 
or more 5^. All the form factors that multiply tensors containing one or 
more factors of 5^, have formulas that do not suffer from the problem of a 
vanishing denominator. These form factors are C 24 , C 35 , C 36 , C 46 , C47, C 48 , 



and C49. The formula for the reduction of C 2 4 is given above in eq. ([)!]). The 
formulas for C35, C36, C46, etc. are given in Appendix B. The rest of the form 
factors, those that multiply tensors without any 5^s in them, are given by 
formulas similar to ([52]) but involving larger matrices. The form factors C 2 i? 
C 22 , and C 23 are given by a set of three simultaneous linear equations, C 3 i 
to C34 by a set of four equations, and C41 to C 45 by a set of five equations. 
The relevant formulas appear in Appendix C. 
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The reduction formulas for these form factors have denominators that are 
proportional to powers of Y { . As discussed earlier, when only one of the Yj's 
is equal to zero no problem arises but when both Y% and Y 2 are equal to zero, 
the limit of these rational functions has to be obtained using l'Hopital's rule. 
We discuss two separate cases here: 

• At least one of the Xij is not zero : In this case the limit is easily 
obtained using l'Hopital's rule or Taylor series expansion where the 
differentiation is performed with respect to one of the mass arguments. 

• All of the Xij are zero : This condition arises if and only if all three of the 
external momenta squared are equal to zero. In this case the Yi cannot 
be shifted away from zero by shifting one of the mass arguments. The 
reduction formulas for this case are obtained as follows. The derivation 
for Co is given below. The rest of the form factors are obtained in the 
same manner. To obtain the reduction for Co, we go back to eq.(|38|). 
When pi = p| = p\ = 0, note that every element of "X -1 " is of the 
0/0 form. Applying l'Hopital's rule on every element of AT -1 , using p\, 
p\, and p\ as independent variables and differentiating with respect to 
p\-i eq.(p8|) becomes, 



C n \ -1 ( | > \ ( R 1 



C 12 / PvP2 \ \ J \ R 2 



(53) 



Since p\-p 2 = and G\\ and Ci 2 are assumed to exist we must have 
that Ri = and R 2 = 0. Substituting for R 1 and R 2 one obtains, 

C = ^{B (2, 3) - B (l, 3)} = Ub (1, 3) - B (l, 2)} (54) 
h h 

Substituting pi = p\ = pi = in (|54] ) one obtains, 

Co(0 M ml,mlml) = ^(O^^mg) - B (0,mf,m§) (55) 



mj — mr, 



B o (0, ml, ml) — B (0, m\, m\ 

2 2 
m 2 ~ m 3 



(56) 



If in the above formulas the denominator vanishes, the limit is obtained 
by differentiating with respect to one of the masses. 
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6.1 Example 

The 1-loop QED contribution to the anomalous magnetic moment (g — 2) 
of the electron was first calculated by Schwinger p0 |. Because of the many 



singularities involved this is one of the most pathological cases that can be 
encountered. The result can be derived in terms of 3-point form factors to 
be 



g — 2 /mev 
2 ~ ~ V 2tt) 



(C u + C l2 ) (57) 



in which the arguments of the form factors are Cij(—m 2 , 0, — m 2 ; A 2 , m 2 , m 2 ) 
where m is the electron mass and e is its charge. Although absent in the 
final result, infrared divergences may appear in intermediate steps and are 
regulated by giving the photon a mass, A. The given set of momenta leads to 
a vanishing Gram determinant and for the given masses causes the extended 
method, described in ref. ||, to break down. In the method described here 
the denominators Y,i in eq.s (fH ) and (f45|) both vanish. To alleviate this 



m 3 can be shifted away from m and the reduction performed in the usual 
way. l'Hopital's or Taylor series expansion is then used to take the limit 
m 3 — > m. This same approach works when applied to the method of ref. || 
and involves differentiation with respect to the mass arguments of B by 
means of the formulas given in Appendix A. The photon mass, A can then 
be safely set to zero. Note that the form factors Cu, C22 and C23 all blow 
up as 7713 — > m, even with a finite photon mass, but that none of these occur 
in the expression for the physical quantity (g — 2). Following this procedure 
one obtains 

C n + C12 = ^2 i 1 + ™ 2 ' m ') ~ fi o(-m 2 ; 0, m 2 )] (58) 
using the result 

B (-m 2 ;0,m 2 ) = B o (0; m 2 , m 2 ) + 2 (59) 

this gives 

0-2 = a_ 
2 ~ 2tt 

in agreement with the well-known result. 
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7 Reduction of D XJ Form Factors 

The derivation of the reduction formulas for the D^ form factors when the 
kinematic determinant, which we again call T>, vanishes, is obtained in an 
identical manner as for the CV, form factor reductions. The kinematic deter- 
minant, T>, is given by, 

V = pIp 2 2 p 2 3 - pl(p 2 -p 3 ) 2 - pl(p 1 -p 3 ) 2 - pl(pi-p 2 ) 2 + 2(p v p 2 )(p 1 -p 3 )(p 2 -p 3 ) 

To obtain the formula for D Q when T> = 0, we start with the general formula 
for the reduction of the Dij form factors. 

Xn X12 X13 

x-2i x 22 x 23 | ( R 2 I (60) 

^31 ^32 ^33 






(Ri\ 


)l 


R 2 




KRs 



where, 

Ri = {/iDo + Co(l,3,4)-C (2,3,4)}/2 (61) 

R2 = {/ 2J Do + C , o(l,2,4)-C , o(l,3,4)}/2 (62) 

#3 = {/3A> + C (l,2,3)-C (l,2,4)}/2 (63) 

p 2 22 r 2 2,22 r 2 22,2 

/i — % Tfi 2 — Pi i f 2 = m 2 - m 3 + p 1 - p 5 , J3 = m 3 -m 4 -p i + p 5 
The elements of the co-factor matrix, X, are, 

XU X 12 Xjs \ I p\p\ - p 2 23 Pi 3 P23 - Pl2P 2 3 Pl2P23 ~ PlPl3 
X21 X 22 X 23 = P13P23 ~ P\2P% V\P\ ~ P^ Pl2Pl3 ~ p\P23 
X 3 1 X 32 X33 J \ pi 2 p 23 - p\pi3 P12P13 ~ P 2 lP23 P\P\ - P 2 12 

(64) 

where pij is used to represent pi ■ pj. Using the same argument as in section 
H, when V = 0, in order for Du, Di 2 , and Di 3 to exist, the following relation 
must hold. 

D = ^{-x i3 C Q (l, 2, 3) + (x i3 - x i2 )C (l, 2, 4) 



+ (x l2 - x iX )C Q {l, 3, 4) + XiiC (2, 3, 4)} (65) 
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where, 



Yi = xnft + x i2 f 2 + x i3 f 3 



and i can take on the values 1, 2 or 3. 

The reduction formulas for the and higher form factors when D = 
are obtained as in the case of the Cy form factors. Here again the In- 
form factors that multiply tensors containing one or more S^s, are given 
by formulas which do not suffer from a vanishing denominator for any set 
of momenta or masses. The rest of the form factors are obtained by solving 
simultaneous linear equations. These form factors have denominators which 
are proportional to powers of Yj. 

As in the case of the Cy form factors, when all of the Y, are zero, the 
limit of the rational functions are obtained using l'Hopital's rule or Taylor 
series expansion. Again there are two cases to consider: 

• At least one of the are non-zero: The limit is obtained using l'Hopital's 
rule where the differentiation is performed with respect to one of the 
mass arguments. 

• All the Xij are zero: This case is treated as for the form factors, 
but there is one difference to note. Recall that in the case of the CV, 
form factors, all the vanished if and only if all the external momenta 
squared were zero, but here all the x^ could vanish even when none of 
the external momenta squared are zero. An example of this is when all 
the external momenta are proportional to one another. Again we use 
the same technique described in the earlier section. Apply l'Hopital's 
rule to each element of "X -1 ". Using pf, p\, p|, p 12 , P13, and p 23 as 
independent variables and differentiating with respect to p\ one obtains, 



If the relation ( |66"1) is again of the 0/0 form, ie. if p\ = p\ = p 23 = 0, 
then the differentiation is performed again, but this time with respect 




(66) 



to, say pi, to obtain, 




(67) 



21 



Proceeding with the assumption that flGTf ) is not of the 0/0 form, eq.flBTiD 
becomes, 




Vlvl - Pis 





(Ri\ 


)l 


R2 




\Rs 



Vl -P23 I I R2 I (68) 
-P23 



Since all the x„ are zero, p\p\ — pis = an d therefore we must have, 

pjR 2 - P23R3 = and P23R2 - v\R3 = 
Substituting for R 2 and R 3 one obtains the expression for D . If (F6| 



is of the 0/0 form, then using fl67|) we get i? 3 = 0, and hence the 



expression for D . The expression obtained for D will be in terms of 
Co, but these Co's necessarily reduce to -Bo's. That is, the vanishing of 
all the sufficient condition for the reduction of Dq, and hence 

all D form factors, to B^s. 



8 Alternate Derivation of the Reduction For- 
mulas 

The formulas obtained above for the reduction of the and form factors 
when their respective kinematic determinants vanish, may be obtained in 
another manner described below. 

Consider the derivation of the reduction formulas for the CV, form factors. 
The reduction for Cq is obtained as before but for the C\ form factors, the 
derivation proceeds in an identical manner as in section ^| only up to eq. flSOp . 
The C24 form factor can be handled differently. Instead of using eq.([H]), one 
can use the relation, 

C 24 = ^{-x l2 B 22 (l, 2) + (x l2 - x a )B 22 (l, 3) + x tl B 22 (2, 3)} (69) 

Eq.(^) is only valid when T> = 0, and it is derived from the general formulas 
for the reduction of C 35 and C 36 form factors, given by [0], 

( C 35 \ X( \{f\C 2i + B 22 (l, 3) - 5 22 (2, 3)} \ 
V C 36 J V\ l{f 2 C 2i + B 22 (l,2)-B 22 (l,3)} J 
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When ( |S9D is substituted for C 2 4 in eq.s (^) and (|5D|), the formulas for Cn 
and C12 are obtained in terms of By form factors which can be eventually 
reduced to B 's. 

The reduction of the Dij form factors proceeds in the same manner. The 
advantage of this method is that one obtains formulas that are not coupled. 
That is, one does not have to solve simultaneous equations to obtain formulas 
for the individual form factors, as was necessary in the previous method. The 
disadvantage here becomes evident when this method is used for obtaining 
the formulas for the higher form factors. Note that to obtain the formulas 
for the C\j form factors we had to use the general formulas for the reduction 
of C^j form factors. In the formula for the reduction of the CV, form factors, 
C35 and C36 show up. From the general formulas for the form factors 
one obtains a formula for C35 and C36 in terms of Bij form factors and the 
quadratically divergent C49. Then from the general formulas for C&j form 
factors one can obtain a formula for C49 in terms of By form factors alone. 
Specifically the By form factors that appear in the formula for C49 is the 
quartically divergent form factor, _B 43 . That is, the quadratically divergent 
C49 is obtained as the difference of quartically divergent B i3 . Substituting 
this formula for C49 into the formulas for the C35 and C36 form factors, 
gives C35 and C36 purely in terms of By form factors which can then be 
substituted into the formulas for the C 2 j form factors to give C 2 i, C22 and 
C23 in terms of By form factors alone. Similarly, to obtain the formulas 
for the C%j form factors one has to drive general formulas for the Cy, form 
factors, which contain the highly divergent B^. The situation is even worse 
for the C*4j form factors and similar problems are faced in the case of the -De- 
form factors. 

Therefore the previous derivation is simpler since to obtain the formulas 
for Cy form factors, we only have to look at the general formulas of the Cj + i j 
form factors, and likewise for the Dy form factors. Also, one does not have 
to deal with the highly divergent form factors which show up in this alternate 
derivation. But no matter how the reduction formulas are derived, once the 
formulas are reduced to 5 's and C 's, they take on unique forms. 

9 Comparison with Earlier Work 

The method employed in ref.s|3], |J [L0| to treat the reduction of scalar in- 
tegrals when the Gram determinant vanished yielded spurious logarithmic 
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divergences that were observed to cancel in physical results. Thus, for exam- 
ple, the logarithmic divergences that appeared in the expressions for G\\ and 
C12, canceled in practical situations when the C M ... were constructed and the 
indices were contracted with the appropriate Lorentz tensors. In this section 
we will show how the cancellation occurs, list all the situations under which 
the cancellation necessarily occurs and also give the situation under which 
this cancellation might not occur. 

The formula given for the C\j form factors in Ref [[3] is, 

C n = a 12 B 1 {l,2) + a 13 B 1 {l,3)-a 23 B {2,3) (71) 
Ci 2 = a 13 Bi(l, 3) + a 23 £ 1 (2, 3) (72) 

where the 

ai2 = -x 22 /Y 2 , a 13 = (x 22 - x 2 i)/Y 2 , a 23 = x 2 i/Y 2 

Comparing eq.s ( |71~D and (72) with the complete expressions given in eq.s ( |49| ) 
and (|50D , it is evident that the terms involving the logarithmically divergent 
C24 are missing from the expressions in eq.s ([Fl]) and (|72|). 

The reason why the terms involving C24 cancel in practice is shown below. 
If one constructs = pi^Cu +p 2fl Ci 2 using Cn and C\ 2 from eq.s (^) and 
(|50|), then one obtains, 

C M = — ( -^1^2-81 (1,2) + {pin +p 2 ^)(x i2 -xu)Bi(l,3) 

1 i 

+ p 2lJl x il B 1 {2, 3) - p lfM x a B (2, 3) + 2(p ltl x n + p 2fl x i2 )C 2i ) (73) 

Note that the term involving C24 is proportional to the light-like vector 

Pfx = P\»Xix + p 2li x i2 (74) 

We note that for all higher rank tensor integrals, C^ v , C^x, etc., the terms 
that are missing from the expressions in ref. || are always proportional to the 
vector P a , where a is one of the indices of the tensor integral. 

Remembering that at the outset we have assumed V = 0, we note the 
following properties of the vector P M 

• P M = if p\ = ap 2 for some real a. 

• pip-Pp = 0. 
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• P2it-Pft = fl- 



it was stated in section [| that when T> — 0, it is necessarily true that 
Pi = ap2, unless at least two of the three external momenta, were space- 
like. Therefore, as long as the external momenta are not space-like, P M is 
identically zero. If at least two of the external momenta are space-like, P M 
is not identically zero, but if it is dotted with any of the external momenta 
or any linear combination of them, it vanishes. Under either of these condi- 
tions, which is almost always the case, the term involving C24 vanishes. The 
only condition under which this term can survive is if at least two of the 
external momenta are space-like and it is contracted with a vector which 
is linearly independent of the external momenta of the three point function. 
Such conditions may occur for the process e + e~ — > Ue^eJ produced by t- 
channel exchange of of a W^-boson. There the one-loop corrections to the 
W + W~^j vertex can yield a form-factor having two external space-like mo- 
menta (i.e. those of the W's) and one external light-like momentum (i.e. that 
of the external photon). In such the method used in ref.s|| ||, [10| 



might yield incorrect results but would be indicated by the non-cancellation 
of the extra spurious diverges. In the program LERG-I this would also be 
indicated by the appearance of arbitrary constants in physical results. 

In the case of the Dij form factors, the situation is similar. The "missing 
terms" are always proportional to the vector pi^Xn + pip^a. + P^Xis- The 
properties of this vector are identical to the one above. It is identically zero 
if any two of pi^ and p^ are linearly dependent, and it vanishes if it is 
dotted with any of its external momenta, and hence any linear combination 
of them. Here again, the only condition under which the "missing terms" 
can survive, is if all three external momenta, pi^, p2^-, and p^^, are linearly 
independent and if D^... is not contracted with any vector which is linearly 
dependent with the external momenta of the four-point tensor integral. 



10 Summary and Conclusion 

We have extended the method for the reduction of Lorentz tensor form fac- 
tors to scalar integrals, introduced by Passarino and Veltman, to regions of 
parameter space including those of vanishing Gram determinants. No at- 
tempt has been made to handle numerical instabilities that may arise and 
the region of vanishing Gram determinant is approached since this problem 
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as been amply discussed in the literature. If required the methods described 
here can be used in this critical region by employing a Taylor series expansion. 

The technique presented here may also be applied in the reduction of five- 
point and higher functions ]2T| when their Gram determinants vanish. This 



is necessarily the case for six-point and higher functions in four space-time 
dimensions. 

In the foregoing it is implicitly assumed that the calculations are per- 
formed in a covariant gauge. In axial gauge, for example, new denominators 
appear in the Feynman integrals that are not immediately amenable to the 
methods described here. 

Once again it should be emphasized that many of the formulas that have 
been derived in this paper are valid only for the divergent and finite parts 
of the integrals and do not necessarily apply to the terms of 0(n — 4) and 
higher. Thus the results cannot be used for the insertion of subdiagrams to 
form 2-loop or higher diagrams. 

The reduction of tensor form factors to scalar integrals can also be im- 
plemented at the 2-loop level but new complications arise. Whereas at 
1-loop there is just a single topology for each of the 2-, 3- and 4-point scalar 
integrals, this is no longer true at 2-loops. For the 2-point function at 2-loops 
there are 20 possible topologies for the scalar integrals. Some of these can 
be calculated in closed form [^J but it is known that general case is not 
expressible in terms of polylogarithms and so are probably best calculated 
numerically. 



Appendix A 



The 2-point form factors, Bij, are defined following Passarino and Veltman 
fTTTl via the relations 



Bq\ B^ B^; B^; B tiuaf3 (p; mi, m 2 ) 



d n q 1; q^; q^q u ] q^q v q a ] qy.q v q a qp 
in 2 [q 2 + m 2 ][(q + p) 2 + m|] 



By invoking Lorentz covariance one may define the form factor, 

B^ = p^Bx 

Here and in what follows the arguments of the B^s are taken to be Bij(p 2 ; m 2 , m\) 
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Similarly, 



B^ a = p^p u p a B zl + {p5}^ a B z2 
Bfj.ua/3 = P/jPvPaPpBa + {pp5} nvapB^ + {55} ^pB^ 

The braces are used to denote the symmetrized product in the Lorentz in- 
dices. The form factors Bn (we take Bq = B$i, B\ = Bn) are all logarith- 
mically divergent as signaled by their having poles at even integer n > 4. 

A and B i2 are quadratically divergent with poles at even integer n > 2. 
i?4 3 is quartically divergent with poles at even integer n > 0. 

Expressions for the derivative of Bq with respect to its momentum argu- 
ment can be found in ref.JJ in the form ([LT]) and various special cases and 
recursion relations are to be found in ref. ||. Here we give the expressions 
for the derivatives with respect to both the momentum and mass arguments. 
Neither derivative exists at the physical threshold, p 2 = —(mi + m 2 ) 2 . 

Away from threshold 

dB . 2 2 2n P 2m i + P 2m l +mf - 2m\m\ + m\ , 2 2 2 
— (p*;ml,n4) = ^ B (p ;m 1 ,m 2 ) 

2 {p 2 -m\ + mf) 2 2 

+ m i 2~F> B » °i m i> m i 

p A D 

+m 2 2 (p2 ~f n +mi) ^ (0; ml m\) (A.l) 
p z D 

(p 2 — m\ + ml) (p 2 — ml + m\) 
p^D 

provided the denominator p 2 D ^ 0. Here 

D = (p 4 + m\ + m\ + 2p 2 m\ + 2p 2 m\ - 2m\m 2 2 ) (A.2) 
= [j) 2 + (mi+m 2 ) 2 ][p 2 + (mi-m 2 ) 2 ]. (A.3) 



For p 2 = and m 2 7^ m 2 

^0/ 2 2 2\ 

(P ;m 1 ,m 2 ) 



dp 



m\ + m 2 



p 2=o 2 ( m I - m 



1\2 



2 2 

'" I; " :> [B (0; m\,m\)- B (0; m 2 ., m 2 2 )] (A.4) 



(m\ — m 2 ) 3 
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which reduces to 



9B 2 2 2 , 
— (p ;m 1 ,m 2/ 



p 2 =0 



6m 2 



(A.5) 



when m 2 = m 2 = 



For the special case, p 2 = — (mi — m 2 ) 2 , for which D vanishes 

. 2 2 2\ 

(p ; m ls m 2 ) 

p 2 =—(m, \ - 

[S o (0;m?,m?)-S (0; 



dp 



^2 ^ 



p 2 =— (mi -ni2) 2 

(m 2 — rrij) 



4^2 



2 2 



p 2 =—(rrii—m2) 2 

(A.6) 



Provided D/0 the derivative with respect to m\ is 



2 2 J , P ^1 m 2 D / 2 2 2\ 

(p ;m 1 ,m 2 ) = 5 (p;m 1 ,m 2 ) 



/ 2 2 2 



p 2 + m\ + m\ 2 2\ , 2m 2 n / n 2 2N 
7j A)(0; mj, mj) + —B (0; m 2 , m 2 ) 



{p 2 + m\- ml) 



(A.7) 



When p 2 = —(mi — m 2 ) 2 

/ 2 2 2 

-(p ;m 1 ,m 2 



/ 2 2 2^ 



p 2 =— (mi— m2) 2 
1 



mi (mi — m 2 ) 2p 2 



[B (0; m 2 , m 2 ) - £ (0; m|, m 2 . 



p 2 =— (mi— rri2) 2 

(A.8) 



For p = and mf 7^ m| eq. ( A.7 ) becomes 



^0 / 2 2 2^ 

^ ! (p;m 1 ,m 2/ 
1 



p 2 =— (mi— m2) 2 



(m 2 — m 2 ) 



1 2 2 2^2 I 5 o(0; m 2 , m 2 ) - £? (0; m 2 , m 2 )] (A.9) 
(wij — m\Y L J 



28 



which reduces to 



9B , 2 2 2 \ 
— (p ;m 1 ,m 2 ) 



p 2 =0 



2m 2 



(A.10) 



for m\ = m 2 = m 2 . 

The partial derivative with respect to m 2 , can be obtained from the fact 
that B is symmetric in its mass arguments. 



Appendix B 



The following are the reduction formulas for the CV, form factors that multi- 
ply tensors which have one 5^ u in them. These form factors are logarithmi- 
cally divergent. 



C24 = 

C35 — 

C36 — 

C46 = 

C47 = 

C*48 — 



~\ + \{hC u + f 2 C 12 + 2m\C Q - B (2, 3)} 
I + \{fiC 21 + f 2 C 23 + 2mlCn + B (2, 3)} 



1 1 

6 + 2 



{/iC 23 + f 2 C 22 + 2m\C 12 - ^(2, 3)} 



-\ + \{hCsi + f 2 C 33 + 2m\C 2l - B (2, 3)} 
~ + \{f\Cu + f 2 C 32 + 2m\C 22 - B 21 (2, 3)} 
~ + l{fiC 33 + f 2 C M + 2m\C 23 + B,(2, 3)} 



All the logarithmically divergent CV,- form factors given above can be written 
using a single formula, with a change in notation, as follows, 



n 



where, 



m + n 



[C mn l\ + -{flC{m+l}nO + f 2 C m {n+l}0 



+2m\C mnQ + (-l) m+1 J B n0 (2,3)} 

( _]\m+n+l 

[C mnl ] = lim(ra - 4)C mn i = ■ — — — — 

n->4 (m + n + 2){n + 1) 
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This notation will become evident with an example. 

Cm, = PlnPluC 2 00 + P2pP2uCo20 + {PlP2}nvC lw + tf^Cooi 

The first index of C counts the number of pi's in the tensor that it multiplies, 
the second index, the number of p 2 s and the third index, the number of <5's. 
In the case of the Bij form factors, the first index counts the number of p's, 
and the second index, the number of 5's. 

Finally the only quadratically divergent form factor C49 is given by, 



C 



19 



\fiC 35 + l -f 2 C m + m\C 24 - B 22 (2, 3) + ^p 2 2 B 21 (2, 3) 
+ \(pI -m 2 2 + ml)B 1 (2, 3) - ±m 2 2 B (2, 3) 
+^(P? - Pi + Pi + 4 ™i ~ 2m 2 2 ~ 2™%) 



The Dij form factors that multiply tensors that have one 5^ in them are 
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given by, 

Am 
-D312 
-D313 
-D416 
-D417 
-D418 
-D419 
-D420 
-D421 



= -1 
1 

~2 
1 

~2 
1 

~2 
1 

~3 
1 

~3 
1 

~3 
1 

~3 
1 

~3 
1 

~3 



^{fiDu + f 2 D 12 + f 3 D l3 + 2m\D Q - C (2, 3, 4)} 
^{/i^2i + / 2 ^24 + / 3 ^25 + 2mjD u + C (2, 3, 4)} 



2{A^24 + / 2 ^22 + / 3 ^26 + 2m^ 12 - C n (2, 3, 4)} 



{fiD 25 + / 2 £>26 + / 3 ^23 + 2m?L> 13 - C 12 (2, 3, 4)} 
UiD 31 + / 2J D 34 + f 3 D 35 + 2m\D 21 - C (2, 3, 4)} 
{/i£>36 + /2A2 + f 3 D 38 + 2m\D 22 - C 21 (2, 3, 4)} 
{A^37 + /2A59 + hD 33 + 2m?L> 23 - C 22 (2, 3, 4)} 
{/i^ + / 2 A»6 + f 3 D 310 + 2mjD 24 + C n (2, 3, 4)} 
{fiD 35 + f 2 D 3W + / 3J D 37 + 2m\D 2b + C 12 (2, 3, 4)} 



^{/lAno + / 2 ^38 + / 3 ^39 + 2m\D 2% - C 23 (2, 3, 4)} 



This set of form factors given above can be written using a single formula, 
with a change in notation, as, 



D 



mnul 



m + n + u + 1 



2 {/l-^{m+l}nuO + /2-D m {n+l}«0 + /s^mnfu+lJO 

+2m?£> mnri0 + (-l) m+1 CU,(2,3,4)} 



The notation is as described in the case of the C form factors. The first, 
second, and third index of the form factor count the number of pi's, 
p 2 's, and p 3 's respectively, in the tensor that it multiplies, and the last index 
counts the number of 5's. 

Finally, the only logarithmically divergent form factor, -D4 22 , is given 
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+2{ Pl -p 3 )D 43a + 2( PrP3 )D 421 + m\D 21 - C 24 (2, 3, 4)} 



Appendix C 

The following are the reduction formulas for the CV, form factors when T> = 0. 
The formulas for Co and the C\j form factors are given in section |6|. The 
formulas for the Cy form factors that multiply tensors that have one or more 
5^ are given in Appendix B. The rest of the Cy form factors are given here. 

C 2 j Form Factors 

/ 2x a h + 3Y t 2x n f 2 \ / C 21 \ /R_2i\ 

xah 4F 4 x a f 2 \ [ C 23 ) = \ R 23 ) 

\ 2x i2 ft 3Yi + 2x l2 f 2 j \C 22 j \ R 22 j 

R 21 = -AmlxnCu - 3x i2 B 21 (l, 2) + 3(x i2 - Xn)B 21 (l, 3) 

+XaB (2,3) - -xn 

R 23 = -2m\{x i2 C n + x a C 12 ) + 3(x i2 - xn)B 21 (l,3) - 2xnB u (2, 3) 

-x i2 B (2, 3) - xn/3 - 2x i2 /3 

R 22 = -4mlx i2 C 12 + 3(xi2- Xa)B 21 (l,3) + 3x a B 21 (2,3) 

2 

+2x i2 B u (2,3) - -x i2 
Csj Form Factors 

/ 3x a fi + 4Yi 3x a f 2 \ 

^i 2 fi Xnfi + 5Yi 2xi\f 2 „ _ =- 








2x i2 fx 




5Y; + ^^2/2 

3x i2 fi 



Xnf 2 

m + 3x i2 f 2 J 



x C 3 = R 3 



V 
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where, 

{C3} T = ( C31 C33 C34 C32 ) {R3} T = ( i?3i i?33 i?34 i?32 ) 

i2 3 i = -6m^ a C 2 i- 4^2^31(1, 2) +4(^ 2 -^1)^31(1, 3) 

/ x 3_ 
-2x a -B (2, 3) + -xn 

R 33 = -2m\(x i2 C 21 + 2 Xil C 23 ) + A(x i2 - xa)B 31 (l, 3) + 2x il B 1 (2, 3) 

. . 1_ 1_ 
+x i2 -B (2, 3) + -x a + -x i2 

E 34 = -2m\{2x i2 C 23 + XiiCaa) + 4(x i2 - ^^(l, 3) - 3x il B 21 (2, 3) 

-2x i2 B x {2, 3) + -x a + -xe 
2 

#32 = -6m?C22 + 4(x i 2-x i i)S3i(l,3)+4x i iS 3 i(2,3) 
+3^2^21(2,3) + -x i2 
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C 4 j Form Factors 



2x iX h + 6Yi 





3xnf 2 



\ 













7Y t 2x iX j 2 
4x i2 /! 






xaf 2 



V 



+ 4x i2 / 2 / 

X C4 = R4 



where, 



{C 4 } 
{R 4 } 



T 




T 



i?4i = -8mlx il C 31 -5x i2 B il (l,2)+5(x i2 -x a )B A1 (l,3) 

g 

+Xii5 (2,3) - -Xn 
5 

i?43 = -2m\{x i2 C 3X + 3xiiC 3 3) + 5(x i2 - x^B^l, 3) - 2x il B 1 (2, 3) 

/ x 3_ 2_ 
-x i2 B (2,3) - -xa - -x i2 
5 5 

-R45 = -4mi(^ 2 C 33 + XnCte) + 5(x i2 - x^B^l, 3) + 3^1-821 (2, 3) 

+2^2^(2, 3) (2^i + 3^ 2 ) 
15 

i? 44 = -2mi(3xi 2 C 34 + x a C 32 ) + 5(x i2 - x a )S 41 (l,3) - 4S a S 31 (2,3) 

. . 1 _ 2_ 
-3a;j 2J B 21 (2,3) - — x a - -x i2 
10 5 

-R42 = -8miX i2 C 32 + 5(x i2 - ^1)541(1,3) + 5xii54i(2, 3) 

2_ 

+4x i2J B 3 i(2,3) - -x i2 
5 

Appendix D 

The following are the formulas for the form factor reductions when V = 0. 
The formula for the reduction of Do is given in section |7|. The formulas for 
the reduction of all Dij form factors which multiply tensors that have one or 
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more in them are given in Appendix B. The formulas for the rest of the 
form factor reductions are given below. 

Dij Form Factors 

Yi + xnfr x n f 2 xnf 3 \ / D u \ / R n \ 

X i2 fl Y, t + X i2 f 2 x i2 f 3 ^12 = -R12 

xah x i3 f 2 Yi + x i3 f 3 J \ D 13 J \ R 13 J 

R n = -2m 2 1 x il D -x i3 Cn(l,2,3) + (x i3 -x t2 )Cn(l,2,A) 

+ (x l2 -x il )C u (l,3A) 
R 12 = -2mlx i2 D -x i3 C 12 (l,2,3) + (x i3 -x i2 )C 12 (l,2,A) 

+ (x t2 - x^Cnil, 3, 4) + x tl C n (2, 3, 4) + x i2 C (2, 3, 4) 
R 13 = -2mlx i3 D + (x i3 -x i2 )C 12 (l,2,4:) + (x i2 -x il )C 12 (l,3,4:) 

+x ll C 12 (2,3,4)+^ 3 C (2,3,4) 

The D 2 j and higher D^ form factor formulas involve large matrices which 
do not fit neatly in the page and so the equations are not written in the matrix 
form. 

D 2 j Form Factors 

(Yi + Xnfi)D 21 + x a f 2 D 24 + x a f 3 D 25 = 

- 2m\x lX D xx - x i3 C 21 (l, 2, 3) + (x i3 - x i2 )C 21 (l, 2, 4) 

+ (x i2 -x il )C 21 (l,3,A) 



x i2 fiD 21 + (3Yi - x i3 f 3 )D 2A + x i2 f 3 D 25 + x a f 2 D 22 + x a f 3 D 26 = 

- 2m\{x i2 D 11 + x n D 12 ) + 2(x l2 - x a )C 21 (l : 3, 4) - 2x i3 C 23 {\, 2, 3) 

+ 2(x i3 - x i2 )C 23 {\, 2, 4) - x il C 11 (2, 3, 4) - x i2 C (2, 3, 4) 



x i3 fiD 21 + x i3 f 2 D 24 + (3Yi - x i2 f 2 )D 25 + Xnf 2 D 26 + x a f 3 D 23 = 
- 2ml{x i3 D ll + x a D 13 ) + 2{x l3 - x i2 )C 23 (l, 2, 4) 

+ 2(x l2 - x n )C 23 {\, 3, 4) - x a C 12 (2, 3, 4) - x l3 C (2, 3, 4) 
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x i2 fiD 24 + (Yi + x i2 f 2 )D 22 + x i2 f 3 D 26 = 

- 2m\x i2 D 12 + (x i2 - x a )C 21 (l, 3, 4) + x a C 21 (2, 3, 4) 

- x l3 C 22 (l, 2, 3) + (x i3 - x i2 )C 22 (l, 2, 4) + x i2 C u (2, 3, 4) 

x i3 fiD 24 + x i2 fiD 25 + x i3 f 2 D 22 + (3Yi - x il f 1 )D 2 6 + x i2 f 3 D 23 = 
- 2m\(x i3 D 12 + x i2 D 13 ) + 2(x i3 - x i2 )C 22 (l, 2, 4) + 2(x i2 - x a ) 

C 23 (l, 3, 4) + 2x n C 23 (2, 3, 4) + x i3 C n (2, 3, 4) + x i2 C 12 (2, 3, 4) 

x i3 fiD 25 + x i3 f 2 D 26 + (Yi + x i3 f 3 )D 23 = 

- 2m\x i3 D 13 + (x i3 - x i2 )C 22 (l, 2, 4) + (x i2 - x a )C 22 (l, 3, 4) 

+ xaC 22 (2,3,4)+x i3 C 12 (2,3,4) 

D 3j Form Factors 

(Yi + x a /i)D 3 i + x a f 2 D 34 + x a f 3 D 35 = 

- 2m\x iX D 21 - x a C 31 (l, 2, 3) + (x i3 - x i2 )C 31 (l, 2, 4) 

+ (^ 2 -x a )C 3 i(l,3,4) 

x i2 fiD 3 i + (3Yi + 2x a fi + x i2 f 2 )D 34 

+ x i2 f 3 D 35 + 2xnf 2 D 36 + 2x a f 3 D 3W = 

- 2m\(x i2 D 2l + 2x n D 24 ) - 3x i3 C 33 (l, 2, 3) 

+ 3(x« - ^ 2 )C 33 (1, 2, 4) + 3(x j2 - x a )C 31 (l, 3, 4) 

+ x, 1 C 11 (2,3,4) + ;r j2 Co(2,3,4) 

x i3 fiD 31 + x i3 f 2 D 34 

+ (3Y + 2x a h + x i3 f 3 )D 35 + 2x n f 2 D 310 + 2x a f 3 D 37 = 

- 2m\(x i3 D 2l + 2x a D 25 ) + 3(x i3 - x i2 )C 33 (l, 2, 4) 

+ 3(x i2 - x a )C 33 (l, 3, 4) + x a C 12 (2, 3, 4) + x l3 C (2, 3, 4) 

2x i2 f 1 D 34 + (4Yi + x i2 f 2 - x i3 f 3 )D 36 

+ 2x i2 f 3 D 3W + x a f 2 D 32 + x a f 3 D 3S = 

- 2mj(2x l2 D 24 + x a D 22 ) - 3x l3 C u (l, 2, 3) 

+ 3(^3 - x i2 )C 34 (l, 2, 4) + 3(:r l2 - xn)C 31 (l, 3, 4) 

-2x a C 21 (2,3,4)-2x i2 Cn(2,3,4) 
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x i3 fiD 3A + x i2 fiD 35 + x i3 f 2 D 36 

+ 4YiD 310 + x i2 f 3 D 37 + x n f 2 D 38 + xnf 3 D 39 = 

- 2m\(x i3 D 2i + x i2 D 25 + x a D 26 ) + 3{x i2 - x a )C 33 (l, 3, 4) 

+ 3fe - x i2 )C u (l, 2, 4) - 2x tl C 23 (2, 3, 4) 

-s«Cii(2,3,4)-s i2 C'i 2 (2,3,4) 

2x i3 fiD 3b + 2x i3 f 2 D 310 

+ (4Yi - x i2 f 2 + x i3 f 3 )D 37 + xnf 2 D 39 + xnf 3 D 33 = 

- 2ml(2x t3 D 25 + x a D 23 ) + 3(x i3 - x t2 )C u (l, 2, 4) 
+ 3(^ 2 - xn)Cu(l, 3, 4) - 2x n C 22 (2, 3, 4) - 2x l3 C 12 (2, 3, 4) 

x i2 fiD 36 + (Y,i + x i2 f 2 )D 32 + x i2 f 3 D 38 = 

- 2m\x i2 D 22 + (x i2 - x a )C 31 {l, 3, 4) + x a C 31 (2, 3, 4) 

- x i3 C 32 (l, 2, 3) + (x l3 - x i2 )C 32 (l, 2, 4) + x l2 C 21 (2, 3, 4) 

Xi 3 fiD 36 + 2x i2 fiD 31{) + x i3 f 2 D 32 

+ (SYi + 2x t2 f 2 + x i3 f 3 )D 38 + 2x t2 f 3 D 39 = 

- 2m\(x i3 D 22 + 2x i2 D 26 ) + 3(x i3 - x i2 )C 32 (l, 2, 4) 

+ 3(x i2 - x tl )C 33 (l, 3, 4) + 3x tl C 33 (2, 3, 4) 

+ x i 3C 2 i(2,3,4) + 2x i2 C 23 (2,3,4) 

2x i3 f 1 D 3W + x i2 f\D 37 + 2x i3 f 2 D 38 

+ (3Yi + x i2 / 2 + 2x i3 f 3 )D 39 + x i2 f 3 D 33 = 

- 2m\(2x l3 D 2& + ^ 2j D 23 ) + 3(^3 - x i2 )C 32 (l, 2, 4) 

+ 3(x l2 - x n )C M {\, 3, 4) + 3x a C 34 (2, 3, 4) 

+ x i2 C 22 (2,3,4) + 2x i3 C 2 3(2,3,4) 

x i3 fiD 37 + x i3 f 2 D 39 + (Fi + x i3 f 3 )D 33 = 

- 2m\x l3 D 23 + (x i3 - x i2 )C 32 (l, 2, 4) + (x i2 - xn)C 32 (l, 3, 4) 

+ x a C 32 (2,3,4)+x, 3 C 22 (2,3,4) 



37 



D 4 j Form Factors 



(it + xnf^D^ + x a f 2 D 44 + x a f 3 D 45 = 

- 2m\x tl D 3X - x i3 C 41 (l, 2, 3) + (x i3 - x i2 )C 41 (l, 2, 4) 

+ (x i2 - x a )C 4 i(l,3,4) 

x i2 fiD 41 + (5Yi + 2x a fi - x i3 f 3 )D 44 

+ x i2 f 3 D 45 + 3x a f 2 D 41Q + 3x a f 3 D 413 = 

- 2m\(x i2 D 31 + 3x a D 34 ) - Ax i3 C 43 (l, 2, 3) 

+ 4(x i3 - x i2 )C 43 (l, 2, 4) + A{x l2 - x a )C 41 (l, 3, 4) 

- x il C 11 (2,3,4)-^ 2 C (2,3,4) 

Xi 3 fiD 4 \ + X{ 3 f 2 D 44 

+ (5Yi + 2x il f 1 - x i2 f 2 )D 45 + 3xnf 2 D 413 + 3x a f 3 D 411 = 
- 2m\{x i3 D 31 + 3x n D 35 ) + A(x i3 - x i2 )C 43 (l, 2, 4) 

+ 4(x t2 - z a )C 43 (l, 3, 4) - x n C 12 (2, 3, 4) - x t3 C (2, 3, 4) 

Xi 2 f\D 44 + (3Yi — x i3 f 3 )D 4W + Xi 2 f 3 D 4 i 3 + Xnf 2 D 4 Q + Xnf 3 D 4 i 4 = 
- 2m\{x l2 D 34 + x a D 36 ) - 2x i3 C 45 (l, 2, 3) + 2(x l3 - x i2 )C 45 (l, 2, 4) 

+ 2(x t2 - x u )C 41 (l, 3, 4) + x n C 21 (2, 3, 4) + x i2 C n (2, 3, 4) 

x% 3 f\D 44 + Xi 2 fiD 4 § + Xi 3 f 2 D 4 iQ 

+ (5Yi + xnfi)D 4 i 3 + x i2 f 3 D 411 + 2x a f 2 D 414 + 2x a f 3 D 415 = 

- 2m\(x i3 D 34 + x i2 D 35 + 2x a D 310 ) 

+ A(x i3 - x i2 )C 45 (l, 2, 4) + A(x i2 - xa)C 43 (l, 3, 4) 
+ 2x n C 23 (2, 3, 4) + x i3 C u (2, 3, 4) + x i2 C 12 (2, 3, 4) 

Xi 3 f\D 4 ^ + X{ 3 f 2 D 4 i 3 

+ (3Yi - x i2 f 2 )D 411 + x a f 2 D 415 + Xnf 3 D 48 = 

- 2m\{x i3 D 3b + x n D 37 ) + 2(x i3 - x i2 )C 45 (l, 2, 4) 

+ 2(x l2 - x n )C 4b {l, 3, 4) + x a C 22 (2, 3, 4) + x i3 C 12 (2, 3, 4) 
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3x i2 /i-D4io + (5Yi + 2x i2 f 2 - x i3 f 3 )D m 

+ 3x i2 f 3 D 414 + x n f 2 D 42 + x a f 3 D 47 = 

- 2m\{3x i2 D m + xnD 32 ) - 4x i3 C 44 (l, 2, 3) 

+ 4(^3 - x i2 )C u (l, 2, 4) + 4{x l2 - x a )C 41 (l, 3, 4) 

-3xaC 3 i(2,3,4)-3x i2 C 21 (2,3,4) 

Xi 3 fiD 4 i + 2xi 2 f\D 4 i 3 + Xi 3 f 2 D 46 

+ (5Yi + x i2 f 2 )D 414 + 2x i2 f 3 D 415 + x a f 2 D 47 + x a f 3 D 412 = 

- 2m\(x i3 D m + 2x i2 D 3W + x a D 38 ) 

+ A(x l3 - x i2 )C 44 (l, 2, 4) + A(x l2 - Xil )C 43 (l, 3, 4) 

- 3x n C 33 (2, 3, 4) - x i3 C 21 (2, 3, 4) - 2x i2 C 23 (2, 3, 4) 

2xi 3 f\D 4 i 3 + Xi 2 f\D 4 \\ + 2xi 3 f 2 D 4 i 4 

+ (5Yi + x i3 f 3 )D 415 + x i2 f 3 D 48 + xnf 2 D 412 + Xnf 3 D 49 = 

- 2ml(2x i3 D 3W + x i2 D 37 + x a D 39 ) 

+ 4(x« - x i2 )C 44 {l, 2, 4) + A(x l2 - x a )C 45 (l, 3, 4) 

- 3x a C 34 (2, 3, 4) - x i2 C 22 (2, 3, 4) - 2x i3 C 23 (2, 3, 4) 

3Xj3/iD 4 ll + 3x i3 f 2 D 415 

+ (5Yi - x i2 f 2 + 2x i3 f 3 )D 48 + Xnf 2 D 49 + x a f 3 D 43 = 

- 2m\{3x i3 D 37 + xnD 33 ) + 4{x i3 - x i2 )C 44 {l, 2, 4) 
+ 4(x l2 - x n )C 44 {l, 3, 4) - 3z a C 32 (2, 3, 4) - 3x i3 C 22 (2, 3, 4) 

x i2 fiD 46 + (Yi + x i2 f 2 )D 42 + x i2 f 3 D 47 = 

- 2m\x i2 D 32 - x i3 C 42 (l, 2, 3) + (x i3 - x i2 )C 42 (l, 2, 4) 

+ (x t2 - x tl )C 41 (l, 3, 4) + x n C 41 (2, 3, 4) + x i2 C 31 (2, 3, 4) 

Xi 3 f\D 4& + 3x i2 fiD 4 i 4 + x i3 f 2 D 42 

+ (5Yj - x a fi + 2x i2 f 2 )D 47 + 3x i2 f 3 D 412 = 

- 2m\(x i3 D 32 + 3x i2 D 38 ) + A(x i3 - x i2 )C 42 (l, 2, 4) 

+ 4(x t2 - x tl )C 43 (l, 3, 4) + 4x a C 43 (2, 3, 4) 

+ x i 3C , 3i(2,3,4) + 3x i2 C 3 3(2,3,4) 
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^i 3 f\D 414 + x i2 fxD 415 + x i3 f 2 D 47 

+ (3Yi - Xii/i)D 4 i2 + x i2 f 3 D 49 = 

- 2m?(x i3 £>38 + ^2^39) + 2(^3 - ^ 2 )C 42 (1, 2, 4) 

+ 2(x i2 - x u )C 45 (l : 3, 4) + 2x a C 45 (2, 3, 4) 

+ x i3 C33(2,3,4)+x i2 C 3 4(2,3,4) 

3^3/1 -D415 + x i2 fiD 48 + Sx i3 f 2 D 4 i 2 

+ (5Fj - x a fi + 2x i3 f 3 )D 49 + x i2 f 3 D 43 = 

- 2m\{?>x i3 D 39 + x i2 D 33 ) + 4{x i3 - x i2 )C 42 (l, 2, 4) 

+ 4(x l2 - x a )C M (l, 3, 4) + 4x ll C 44 (2, 3, 4) 

+ x l2 C 32 (2,3,4) + 3x,3C3 4 (2,3,4) 

x i3 fiD 48 + x i3 f 2 D 49 + (Yi + x i3 f 3 )D 43 = 

- 2m\x i3 D 33 + (x i3 - x i2 )C 42 (l, 2, 4) + (x i2 - x a )C 42 (l, 3, 4) 

+ x a C 42 (2,3,4)+x i3 C3 2 (2,3,4) 
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